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LONG-TERM  GOALS 

The  long-term  scientific  goals  of  this  research  project  are  to: 

1 .  Understand  and  quantify  the  sources  of  error  in  ocean  models  that  fundamentally  limit 
the  practical  predictability  of  the  coastal  ocean  circulation. 

2.  Use  information  about  model  error  to  improve  ocean  circulation  estimates  obtained 
using  weak  constraint  data  assimilation  methods. 

OBJECTIVES 

The  primary  objective  of  the  proposed  research  is  to  develop  a  weak  constraint,  4-dimensional 
variational  (4D-Var)  data  assimilation  capability  for  the  Regional  Ocean  Modeling  System 
(ROMS)  with  application  to  the  California  Current  System  (CCS).  The  CCS  is  of  considerable 
socio-economic  and  strategic  significance  to  the  United  States,  and  ROMS  CCS  has  transitioned 
to  a  near  real-time  analysis  system  in  support  of  the  U.S.  west  coast  components  of  the 
Integrated  Ocean  Observing  System  (IOOS).  This  project  is  therefore  very  timely  given  the 
limiting  nature  of  model  errors  on  coastal  ocean  prediction. 


1 


Report  Documentation  Page 

Form  Approved 

OMB  No.  0704-0188 

Public  reporting  burden  for  the  collection  of  information  is  estimated  to  average  1  hour  per  response,  including  the  time  for  reviewing  instructions,  searching  existing  data  sources,  gathering  and 
maintaining  the  data  needed,  and  completing  and  reviewing  the  collection  of  information.  Send  comments  regarding  this  burden  estimate  or  any  other  aspect  of  this  collection  of  information, 
including  suggestions  for  reducing  this  burden,  to  Washington  Headquarters  Services,  Directorate  for  Information  Operations  and  Reports,  1215  Jefferson  Davis  Highway,  Suite  1204,  Arlington 

VA  22202-4302.  Respondents  should  be  aware  that  notwithstanding  any  other  provision  of  law,  no  person  shall  be  subject  to  a  penalty  for  failing  to  comply  with  a  collection  of  information  if  it 
does  not  display  a  currently  valid  OMB  control  number. 

1.  REPORT  DATE 

30  SEP  2013  2' REPORT  TYPE 

3.  DATES  COVERED 

00-00-2013  to  00-00-2013 

4.  TITLE  AND  SUBTITLE 

Quantifying  the  Amplitude,  Structure  and  Influence  of  Model  Error 
during  Ocean  Analysis  and  Forecast  Cycles 

5a.  CONTRACT  NUMBER 

5b.  GRANT  NUMBER 

5c.  PROGRAM  ELEMENT  NUMBER 

6.  AUTHOR(S) 

5d.  PROJECT  NUMBER 

5e.  TASK  NUMBER 

5f.  WORK  UNIT  NUMBER 

7.  PERFORMING  ORGANIZATION  NAME(S)  AND  ADDRESS(ES) 

University  of  California, Department  of  Ocean  Sciences, 1156  High 

St, Santa  Cruz, CA, 95064 

8.  PERFORMING  ORGANIZATION 

REPORT  NUMBER 

9.  SPONSORING/MONITORING  AGENCY  NAME(S)  AND  ADDRESS(ES) 

10.  SPONSOR/MONITOR'S  ACRONYM(S) 

11.  SPONSOR/MONITOR'S  REPORT 
NUMBER(S) 

12.  DISTRIBUTION/AVAILABILITY  STATEMENT 

Approved  for  public  release;  distribution  unlimited 

13.  SUPPLEMENTARY  NOTES 

14.  ABSTRACT 

15.  SUBJECT  TERMS 

16.  SECURITY  CLASSIFICATION  OF:  17.  LIMITATION  OF 

_ _ _  ABSTRACT 

18.  NUMBER  19a.  NAME  OF 

OF  PAGES  RESPONSIBLE  PERSON 

a.  REPORT  b.  ABSTRACT  c.  THIS  PAGE  Same  OS 

unclassified  unclassified  unclassified  Report  (SAR) 

10 

Standard  Form  298  (Rev.  8-98) 

Prescribed  by  ANSI  Std  Z39-18 


APPROACH 


The  proposed  research  activities  seek  to  quantify  the  amplitude,  structure  and  influence  of  model 
error  in  ocean  analysis  and  forecast  systems  using  two  approaches.  The  first  uses  the  ROMS  4D- 
Var  data  assimilation  systems  (Moore  et  al,  2010a,b,c)  in  combination  with  Bayesian 
hierarchical  modeling  to  identify  sources  of  model  error  and  their  characteristics.  The  second 
approach  uses  Generalized  Stability  Analysis  (GSA)  (Moore  et  al,  2004)  to  explore  the  fastest 
growing  model  error  structures,  how  these  relate  to  different  ocean  circulation  regimes,  and  the 
implied  bounds  on  model  error  growth  and  predictability. 


The  project  involves  three  major  research  tasks: 

Research  Task  #1: 

(i)  A  sequence  of  strong  constraint  4D-Var  experiments  for  ROMS  CCS  spanning  several  years 
will  be  used  to  identify  geographical  areas  where  model  error  influences  surface  forcing 
adjustments  imposed  by  4D-Var.  (ii)  A  surface  flux  BHM  will  be  used  to  identify  periods  when 
model  error  identified  in  (i)  is  a  significant  factor.  The  basis  for  the  CCS  BHM  will  be  the  BHM 
of  Milliff  et  al  (2011)  for  the  Mediterranean  surface  winds  which  will  be  reconfigured  for  the 
ROMS  CCS  domain  and  expanded  to  include  surface  fluxes  of  heat  and  momentum.  The  BHM 
data  stage  will  use  QuikSCAT  surface  vector  winds,  and  COAMPS  surface  winds  and  fluxes, 
while  the  process  model  stage  will  utilize  the  bulk  surface  flux  model  subcomponent  of  ROMS 
and  COAMPS  standard  10  m  atmospheric  boundary  layer  variables,  (iii)  A  second  sequence  of 
weak  constraint  4D-Var  experiments  using  the  surface  flux  BHM  and  strong  constraint  flux 
increments  to  inform  the  model  error  prior  will  be  used  to  build  and  test  various  different 
hypotheses  about  model  error,  (iv)  A  detailed  analysis  of  the  spatio-temporal  corrections  for 
model  error  from  each  weak  constraint  4D-Var  assimilation  cycle  will  be  performed  to  identify 
the  nature  of  the  model  errors. 

Research  Task  #2: 

A  complete  characterization  of  model  error  during  each  4D-Var  cycle  is  provided  by  the 
posterior  covariance.  ROMS  4D-Var  permits  computation  of  the  leading  eigenvectors  (EOFs)  of 
the  posterior  error  covariance  matrix.  Spatio-temporal  variations  in  the  structures  of  the  leading 
EOFs  of  posterior  error  will  yield  posterior  information  about  the  efficacy  of  the  resulting  4D- 
Var  circulation  estimates.  Therefore  a  detailed  analysis  of  the  leading  EOFs  of  the  posterior 
model  error  covariance  matrix  will  also  be  performed  for  the  assimilation  sequences  performed 
during  Task  #l(i)  for  strong  constraint  and  Task  #1  (iii)  for  weak  constraint. 

Research  Task  #3: 

(i)  The  stochastic  optimals  (SOs)  and  forcing  singular  vectors  (FSVs)  of  any  time  evolving 
circulation  can  be  computed  using  the  ROMS  GSA  took-kit  described  by  Moore  et  al  (2004).  A 
systematic  study  of  the  FSVs  and  SOs  associated  with  model  error  will  be  performed  for  ROMS 
CCS  using  hindcasts  for  the  same  period  as  Task  l(i).  Hindcast  initial  conditions  will  be 
generated  using  weak  constraint  4D-Var  during  Task  1.  Since  the  leading  FSVs  and  SOs  are  the 
most  damaging  model  errors  for  the  hindcast  interval,  they  yield  upper  bounds  on  model  error 
growth,  and  the  loss  of  predictability  of  the  circulation  during  each  hindcast  interval  due  to  the 
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growth  of  model  error,  (ii)  The  projection  of  the  posterior  model  errors  diagnosed  from  the  weak 
constraint  4D-Var  experiments  of  Task  #1  onto  the  FSVs  and  SOs  will  also  be  examined. 

WORK  COMPLETED 

In  addition  to  the  Pis,  other  personnel  involved  in  this  project  include:  Dr.  Polly  Smith,  a  post¬ 
doctoral  scholar  hired  to  work  on  this  project;  and  Mr.  Kevin  Smith,  a  graduate  student 
researcher  hired  to  work  on  this  project.  Polly  Smith  left  the  project  in  June  2012  to  take  up  a 
position  in  the  U.K.,  and  was  replaced  by  another  post-doctoral  scholar  Dr.  Emilie  Neveu. 
Unfortunately,  Dr.  Neveu  also  left  the  University  in  early  2013  to  take  up  another  position  in 
France.  In  addition  to  the  loss  of  critical  research  personnel,  the  project  also  ran  short  of  funds 
because  of  the  recent  sequestration  by  the  federal  government.  As  a  result  Mr.  Smith,  the 
graduate  student  funded  to  work  on  this  project,  was  unsupported  during  the  spring  and  summer 
quarters  of  2013.  Despite  these  extenuating  circumstances,  work  has  continued  on  the  project 
during  the  last  fiscal  year,  albeit  at  a  slower  rate.  We  are  actively  looking  to  hire  a  new  post¬ 
doctoral  researcher  to  work  on  this  project. 

Task  #1: 

Work  has  continued  on  Task  1  activities  following  on  from  the  work  led  and  initiated  by  Dr. 
Polly  Smith.  Dr.  Smith  generated  a  sequence  of  strong  constraint  4D-Var  analyses  for  the  period 
2002-2004,  and  in  these  experiments,  the  ROMS  4D-Var  control  vector  comprised  the  initial 
conditions,  surface  forcing  and  open  boundary  conditions.  The  surface  forcing  fields  were 
derived  from  output  from  COAMPS  which  is  known  to  perform  well  in  this  region  (Doyle  et  al, 
2009).  Our  working  hypothesis  is  therefore  that  any  significant  departures  of  the  4D-Var 
adjusted  surface  fluxes  from  the  COAMPS  priors  is  a  strong  indication  of  the  influence  of  model 
error.  The  significance  of  the  4D-Var  surface  flux  adjustments  was  assessed  using  a  surface  flux 
BHM  developed  by  PI  Milliff.  So  far  the  BHM  has  been  developed  only  for  the  surface  wind 
stress,  but  will  be  further  developed  later  to  include  surface  fluxes  of  heat  and  fresh  water.  The 
BHM  provides  an  estimate  of  the  probability  distribution  for  the  surface  winds  which  can  be 
used  to  quantify  the  efficacy  of  the  ROMS  4D-Var  wind  adjustments. 

By  treating  the  COAMPS  surface  flux  fields  as  “perfect,”  a  series  of  model  integrations  spanning 
each  data  assimilation  cycle  and  forced  with  the  4D-Var  adjusted  fields  were  performed  and  used 
to  assess  the  impact  of  model  errors  on  the  resulting  circulation  estimates.  These  calculations 
were  then  used  to  estimate  the  parameters  necessary  for  the  model  error  covariance  matrix  that  is 
required  for  weak  constraint  data  assimilation.  The  parameters  so  derived  were  then  used  for 
Task  l(iii). 

In  order  to  prove  the  concept  of  our  approach,  Dr.  Smith  perfonned  all  of  the  aforementioned 
experiments  using  a  low  resolution  (~30  km  grid-spacing)  version  of  ROMS  CCS  (known  as 
WC13),  which  has  for  several  years  served  as  an  efficient  test  case  for  all  of  the  ROMS  4D-Var- 
based  algorithms.  Unfortunately,  in  this  case  some  significant  issues  relating  to  the  convergence 
of  the  weak  constraint  4D-Var  algorithm  were  discovered  during  several  of  the  4D-Var  cycles, 
calling  into  question  the  validity  of  our  results.  These  issues  appear  to  be  related  to  assimilating 
hydrographic  observations  close  to  the  coast  and  in  relatively  shallow  water.  Despite  repeated 
attempts,  we  were  unable  to  mitigate  the  problem  in  the  low  resolution  WC13,  so  a  decision  was 
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made  to  use  a  higher  resolution  version  of  ROMS  CCS  (known  as  WC12)  with  ~10  km  grid¬ 
spacing.  This  model  has  been  the  workhorse  of  a  recent  30  year  sequence  of  historical  analyses 
for  the  CCS,  and  a  near  real-time  4D-Var  and  prediction  system  currently  run  an  UCSC  in 
support  of  CeNCOOS.  This  model  configuration  is  well  behaved  during  weak  constraint  4D-Var 
experiments,  and  does  not  exhibit  any  of  the  issues  found  in  WC13.  The  sequence  of 
experiments  outlined  above  is  now  being  repeated  using  a  subset  of  years  from  the  30  year 
historical  analysis.  Identification  of  the  nature  of  the  model  errors  described  in  Task  l(iv)  will  be 
explored. 

Task  #2: 

Due  to  the  technical  delays  associated  with  Task  #1,  and  the  lack  of  research  personnel  that 
could  be  supported  on  this  award  during  the  last  fiscal  year,  work  has  yet  to  begin  on  Task  2. 

Task  #3: 

Work  has  continued  on  Task  #3  and  led  by  PI  Moore  in  conjunction  with  a  graduate  student 
Kevin  Smith.  Following  some  pedagogic  applications  of  singular  value  decomposition  (SVD)  to 
a  steady  analytical  zonal  jet  in  a  periodic  channel  using  ROMS,  Smith’s  calculations  have  been 
extended  to  a  complex  time-dependent  case  in  a  larger  domain  with  high  horizontal  resolution. 
The  basic  set-up  of  the  model  is  the  relaxation  of  a  zonally  uniform  meridional  temperature 
gradient  and  the  subsequent  development  of  a  baroclinically  unstable  jet.  The  initial  focus  of 
these  experiments  was  to  explore  the  structure  and  energetics  of  singular  vectors  conditioned  on 
the  growth  of  the  perturbation  energy  norm.  However,  since  the  ultimate  focus  of  the  proposed 
research  is  the  influence  of  model  error  on  ocean  predictability,  a  more  appropriate  choice  of 
norm  for  the  SVD  analyses  is  one  based  on  the  analysis  error  covariance  matrix,  Ea,  of  a  forecast 
initial  condition  and  the  forecast  error  covariance,  E  ,  at  the  forecast  time.  Following  Ehrendorfer 
and  Tribbia  (1995),  (E3)'1  is  used  as  the  norm  at  initial  time,  while  E1  is  used  as  the  norm  at  final 
time.  As  shown  by  Ehrendorfer  and  Tribbia  (1998)  the  appropriately  rescaled,  time-evolved  SVs 
computed  in  this  case  are  the  EOFs  of  the  expected  forecast  error  covariance  matrix.  Because  of 
the  intimate  connection  between  Ea  and  the  4D-Var  cost  function  Hessian,  the  SVs  that  evolve 
into  the  EOFs  are  often  also  referred  to  as  Hessian  singular  vectors  (HSVs)  (Barkmeijer  et  al, 
1998).  During  the  last  year  we  have  continued  to  develop  the  mathematical  framework  and 
methodology  necessary  to  compute  the  HSVs  in  ROMS,  and  have  extended  the  same  ideas  to 
stochastic  forcing  representing  model  error  to  yield  what  we  refer  to  as  Hessian  stochastic 
optimals  (HSOs).  To  our  knowledge  the  idea  of  HSOs  is  new.  Another  aspect  of  our  work  which 
is  new  is  that  we  compute  the  HSVs  (and  HSOs)  only  in  the  subspace  that  is  spanned  by  the  4D- 
Var  increments  to  the  prior  circulation  estimate.  This  results  in  significant  computational 
savings,  and  is  in  contrast  to  related  calculations  in  numerical  weather  prediction  where  the 
HSVs  are  typically  identified  in  the  full  state-space  of  the  model.  An  extensive  sequence  of  HSV 
calculations  have  been  performed  for  the  unstable  baroclinic  jet  and  for  the  CCS  by  Smith  and 
Moore,  and  we  are  in  the  process  of  analyzing  them  with  the  view  of  submitting  a  peer-reviewed 
publication  in  late  2013  or  early  2014. 

We  have  also  been  developing  some  new  ideas  in  relation  to  SVD.  The  first  of  these  involves 
using  the  entire  4D-Var  control  vector  (i.e.  initial  condition  error,  surface  forcing  error,  open 
boundary  condition  error,  and  model  error)  as  the  state-vector  for  the  SV  calculations.  The 
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advantage  of  this  is  that  many  of  the  assumptions  that  are  required  to  compute  the  SOs,  such  as 
separable  space-time  correlations  and  errors  that  are  Markovian  in  time,  are  no  longer  necessary 
so  the  resulting  analyses  will  be  more  general.  We  refer  to  these  new  singular  vectors  as  control 
singular  vectors.  The  second  idea  that  we  are  developing  relates  to  SVD  using  Ea  and  E*  to 
define  the  norms.  The  computation  of  both  covariance  matrices  is  non- trivial,  and  in  numerical 
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weather  prediction  it  is  usually  assumed  that  a  good  surrogate  for  E  is  the  energy  nonn.  It 
remains  to  be  seen  if  this  is  true  in  the  ocean.  However,  as  part  of  this  project  we  have  developed 
a  method  for  computing  both  Ea  and  E*  using  the  adjoint  of  the  entire  4D-Var  system  (Moore  et 
ah,  2012).  By  combining  the  adjoint  of  4D-Var  with  the  GSA  tool-kit,  we  will  be  able  to  use  the 
true  expected  forecast  error  covariance  matrix  as  the  appropriate  norm  for  SVD.  A  key  element 
to  advancing  these  activities  has  been  the  development  and  testing  of  the  tangent  linear  and 
adjoint  of  the  Lanczos  formulation  of  the  Restricted  Preconditioned  Conjugate  Gradient  (RPCG) 
algorithm  for  ROMS.  As  reported  last  year,  the  Lanczos  form  of  RPCG  was  implemented  and 
tested  in  ROMS  as  part  of  this  project  and  is  described  in  Giirol  et  al.  (2013). 

Augmented  Restricted  Preconditioned  Conjugate  Gradient  -  Lanczos  Formulation 

One  additional  and  very  important  research  activity  that  has  been  undertaken  as  part  of  this 
project  is  the  implementation  of  a  Lanczos  formulation  of  the  Augmented  Restricted 
Preconditioned  Conjugate  Gradient  (RPCG)  method  described  by  Gratton  and  Tshimanga 
(2009).  The  implementation  of  Lanczos-RPCG  in  ROMS  was  described  in  last  year’s  annual 
report  and  is  documented  by  Giirol  et  al.  (2013)  for  the  case  of  a  single  outer-loop.  In  the  case  of 
multiple  outer-loops,  the  augmented  form  of  the  algorithm  is  required  to  ensure  that  the  resulting 
posterior  increments  lie  in  the  space  spanned  by  the  previous  outer-loops.  The  augmented 
Lanczos-RPCG  algorithm  has  been  successfully  implemented  and  tested  in  ROMS,  and  will  be 
released  to  the  ROMS  community  in  the  near  future.  This  is  an  very  significant  development  for 
ROMS  4D-Var,  and  overcomes  the  poor  convergence  behavior  of  the  dual  4D-Var  algorithms 
reported  by  Moore  et  al  (2011b)  for  ROMS  but  also  widely  acknowledged  in  the  field.  RPCG  is 
formulated  so  as  to  guarantee  the  same  rate  of  convergence  as  the  primal  fonnulation.  What 
makes  this  so  significant  is  that  weak  constraint  4D-Var  is  only  practical  in  the  dual  fonnulation, 
so  the  fast  convergence  enjoyed  by  the  strong  constraint  primal  approach  using  multiple  outer- 
loops  can  now  be  achieved  during  multiple  outer-loops  of  weak  constraint  4D-Var  also.  An 
example  calculation  is  shown  in  the  Results  section. 

RESULTS 

Hessian  Singular  Vectors 

As  noted  above,  the  relevant  norms  for  the  predictability  problem  for  SVD  are  based  on  the 
analysis  error  covariance  matrix  and  forecast  error  covariance  matrix.  Using  the  latter  to  define 
the  final  time  nonn  and  the  inverse  of  the  fonner  to  define  the  initial  nonn  yields  the  so-called 
Hessian  singular  vectors  (HSVs).  The  rescaled  time  evolved  HSVs  are  the  EOFs  of  the  forecast 
enor  covariance  matrix.  Some  example  calculations  are  illustrated  below. 

Figure  1  shows  the  SSH  component  of  the  EOFs  (i.e.  the  time  evolved  HSV)  for  a  6  day  forecast 
of  a  southern  hemisphere  unstable  baroclinic  jet  in  a  zonal  channel  (1000  km  x  2000  km)  for  one 
member  of  a  sequence  of  identical  twin  experiments.  The  forecast  was  initialized  using  an  ocean 
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state  estimate  derived  using  strong  constraint  4D-Var  and  imperfect  simulated  observations  of 
the  upper-ocean  temperature  every  24  hours,  and  a  data  assimilation  window  of  48  hours.  The 
forecast  SSH  at  the  initial  and  final  forecast  times  is  also  shown,  along  with  the  percentage  of  the 
expected  forecast  error  variance  accounted  for  by  each  EOF.  The  first  3  members  of  the  EOF 
spectrum  account  for  almost  ~60%  of  expected  forecast  error  variance. 
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Figure  1 :  The  upper  left  panels  show  the  sea  surface  height  field  at  initial  and  final  forecast  time  for  a 
representative  6  day  forecast  of  a  southern  hemisphere  baroclinically  unstable  jet  in  a  zonal  channel 
(1000  km  x  2000  km).  The  percentage  of  the  expected  forecast  error  variance  explained  by  each  EOF  is 
also  indicated.  The  EOFs  are  the  time  evolved  HSVs  computed  using  the  inverse  of  the  expected 
analysis  error  covariance  as  the  initial  time  norm  and  the  expected  forecast  error  covariance  matrix  as  the 
final  time  norm.  In  order  to  account  for  the  different  units  of  the  forecast  state-vector  elements,  the  final 
time  forecast  error  is  defined  in  terms  of  the  perturbation  energy  of  the  deviations  of  the  forecast  from  the 
truth. 


Figure  2  shows  an  example  calculation  from  ROMS  CCS  (WC12).  In  this  example  actual 
satellite  and  in  situ  hydrographic  observations  were  assimilated  into  ROMS  CCS  using  strong 
constraint  4D-Var  and  a  8  day  assimilation  window  (see  http://oceanmodeling.ucsc.edu).  In  this 
case  the  first  three  members  of  the  EOF  spectrum  account  for  ~70%  of  the  expected  forecast 
error  covariance. 

At  the  present  time  we  are  actively  exploring  the  link  between  the  properties  of  the  EOF 
spectrum  derived  directly  from  the  HSVs  and  the  actual  forecast  skill  of  the  model,  both  in  the 
case  of  the  baroclinically  unstable  jet  and  for  the  CCS.  This  may  provide  a  means  for  predicting 
the  expected  level  of  forecast  error. 
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Figure  2:  The  upper  left  hand  panel  shows  the  8  day  forecast  sea  surface  height  field  initialized  from  an  8 
day  4D-Var  data  assimilation  cycle.  The  remaining  panels  show  the  day  8  sea  surface  height  component 
of  the  leading  members  of  the  EOF  spectrum  computed  from  the  time  evolved  HSVs,  and  in  each  case 
the  percentage  of  the  expected  forecast  error  variance  explained  by  each  EOF  is  indicated. 


Augmented  Restricted  Preconditioned  Conjugate  Gradient  (RPCG) 

As  mentioned  earlier,  an  important  development  in  the  ROMS  weak  constraint  4D-Var  system  in 
the  last  year  has  been  the  implementation  of  the  Lanczos  version  of  augmented  RPCG  which 
guarantees  that  the  primal  and  dual  formulations  of  ROMS  4D-Var  will  converge  at  the  same 
rate  when  using  multiple  outer-loops.  This  is  illustrated  in  Fig.  3  which  shows  time  series  of  the 
4D-Var  cost  function  ,J  from  ROMS  CCS  (WC12)  from  three  experiments.  In  one  experiment, 
the  primal  form  of  the  strong  constraint  4D-Var  algorithm  is  used  in  conjunction  with  5  outer- 
loops  and  15  inner-loops,  and  the  cost  function  convergence  during  each  outer-loop  is  relatively 
fast  and  monotonic.  The  jumps  that  occur  in /at  the  start  of  each  outer-loop  are  primarily  due  to 
two  sources  of  non-linearity  in  the  model  associated  with  (i)  3ld-order  upwind  finite  differencing, 
and  (ii)  GLS  vertical  mixing.  In  a  second  experiment,  the  data  assimilation  was  repeated  using 
the  dual  form  of  strong  constraint  4D-Var  in  conjunction  with  the  Lanczos  version  of  augmented 
RPCG.  Figure  3  shows  that  in  agreement  with  the  theory,  the  primal  and  dual  cost  functions 
track  each  other  perfectly  across  all  inner-  and  outer-loops.  In  the  third  experiment,  dual  weak 
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constraint  4D-Var  was  applied  in  conjunction  with  Lanczos  augmented  RPCG  to  demonstrate 
the  superior  performance  of  the  system  in  this  case.  Recall  that  weak  constraint  is  possible  in 
ROMS  only  using  the  dual  form  of  4D-Var. 
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Figure  3:  A  time  series  of  the  cost  function  J  versus  the  total  number  of  iterations  (i.e.  number  of  inner- 
loops  x  number  of  outer-loops)  for  three  different  ROMS  4D-Var  calculations  in  the  CCS  (1 0  km 
resolution,  42  o-levels)  for  a  representative  7  day  period  29  March  -  4  April,  2003,  each  using  5  outer- 
loops  and  15  inner-loops.  In  one  case  the  strong  constraint  primal  formulation  was  used  (red  dots),  while 
in  a  second  case  the  strong  constraint  dual  formulation  was  used  in  conjunction  with  the  Lanczos 
formulation  of  the  augmented  RPCG  algorithm  (blue  curve).  As  required  the  dual  formulation  using 
augmented  Lanczos-RPCG  yields  identical  iterates  to  the  primal  formulation.  In  the  third  case,  the  weak 
constraint  dual  4D-Var  algorithm  was  applied  in  conjunction  with  augmented  Lanczos  RPCG  (green 
curve). 


IMPACT/APPLICATIONS 

This  project  contributes  significantly  to  the  functionality  and  utility  of  ROMS,  a  widely  used  and 
important  community  model  and  resource.  ROMS  is  unique  in  that  of  all  the  community  ocean 
models  that  are  available,  it  is  the  only  model  that  possesses  such  a  wide  range  of  4D-Var 
algorithms,  analysis  tools,  and  diagnostic  capabilities.  As  part  of  this  project,  we  have  already 
added  to  the  utility  of  ROMS  by  the  addition  of  the  Lanczos  formulation  of  the  augmented 
RPCG  algorithm,  time  correlations  for  model  error  in  weak  constraint  4D-Var,  and  three  new 
drivers  for  computing  Hessian  SVs,  Hessian  Forcing  SVs,  and  Hessian  SOs. 

TRANSITIONS 

The  new  ROMS  utilities  developed  as  part  of  this  project  will  be  freely  available  from  the 
ROMS  web  site  http://myroms.org  and  will  be  actively  used  and  further  developed  by  other 
research  groups  in  the  U.S.  and  elsewhere  as  user  competence  increases. 

RELATED  PROJECTS 

The  work  described  here  is  closely  related  to  the  following  ONR  supported  projects: 

“A  community  Terrain-Following  Ocean  Model  (ROMS)”,  PI  Hernan  Arango,  grant  number 
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NOOO 14-08- 1-0542. 


“Bayesian  Hierarchical  Models  to  Augment  the  Mediterranean  Forecast  System”,  PI  Ralph  Miliff, 
grant  number  N00014-05-C-0198. 

“Understanding  Predictability  of  the  Ocean”,  PI  Brian  Powell,  grant  number  N00014-09-1-0939. 

“Bayesian  Hierarchical  Model  Characterization  of  Model  Error  in  Ocean  Data  Assimilation  and 
Forecasts”,  PI  Ralph  Milliff,  grant  number  N00014-10-C-0354. 

In  addition,  the  development,  implementation  and  testing  of  the  Lanczos  formulation  of  the  augmented 
RPCG  algorithm  in  ROMS  dual  4D-Var  is  work  that  has  been  done  in  close  collaboration  with 
researchers  at  CERFACS  in  Toulouse.  This  project  has  benefited  considerably  from  direct  interactions 
with  not  only  with  CERFACS  scientists,  but  also  scientists  at  ECMWF. 
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